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A new functional-based approach is developed for the stability analysis of linear 
impulsive systems. The new method, which introduces looped-functionals, con- 



■ siders non-monotonic Lyapunov functions and leads to LMIs conditions devoid 



of exponential terms. This allows one to easily formulate dwell-times results, 
for both certain and uncertain systems. It is also shown that this approach 
may be applied to a wider class of impulsive systems than existing methods. 
Some examples, notably on sampled-data systems, illustrate the efficiency of 
the approach. 
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1. Introduction 

t^i 1 Impulsive systems 0, Sll 0, H E 0, E | are an important class of hybrid 



C> systems admitting discontinuities in the state trajectories, that are governed 

by discrete-time maps. They occur in several fields like epidemiology [23L lij]. 
sampled-data and networked control systems fl7| , forestry |25| , power electron- 
ics [HI, harvesting problems [HI, EH, e t c - Among the wide class of impulsive 
dynamical systems, we may distinguish between systems whose impulse-times 
depend on the system state and those for which the impulse-times are external 
to system and only time-dependent. The latter class may be represented in the 
following form 

x{t) = Ax(t), t e K+\I 
x + (t) = Jx(t), tel (1) 



Email addresses: briatcObsse . ethz . ch, corentinObriat . inf o (Corentin Briat), 
Alexandre.SeuretOgipsa-lab.grenoble-inp.fr (Alexandre Seuret) 

URL: http://www.briat.info (Corentin Briat), 
http://www.gipsa-lab.fr/~-ll~alexandre.seuret (Alexandre Seuret) 



Preprint submitted to Elsevier 



July 10, 2012 



where x, xq G R" are the state of the system and the initial condition, respec- 
tively. The system matrices are possibly uncertain, i.e. (A, J) G A x J where 

A:=co{A 1 ,...,A Na }, A t eR nxn , i = l,...,N A 
J:=co{J u ...,J Nj }, Ji G R nxn , i = 1, . . . ,Nj [Z> 

where Na,Nj 6 N\{0} and co{-} is the convex hull. The state trajectory is 
assumed to be left-continuous and to have right-limits at all time. The nota- 
tion x + {t) denotes the right limit of x(s) as s tends to t from the right, i.e. 
x + (t) = lima;(s). The set of impulsive times I :— {tk}keN is a countable set of 

impulse instants tk+i > t^, k G N and we define the inter-impulse distance as 
Tk '■= tk+i — tk- This quantity is also referred to as dwell-time in the literature. 
We also assume that the sequence {ffcjfceN has no accumulation point in order 
to exclude any Zeno behavior. 

Depending on the structure of the matrices A and J, the system may exhibit 
very different behaviors. In particular, notions of minimal and maximal dwell- 
time can be defined for impulsive systems [l2| , similarly as for switched systems 
TSUI- In the case of impulsive systems, these notions refer to system properties 



such that too large or too short dwell-times destabilize the system. In the case 
of periodic impulses with period T > 0, the problem essentially reduces to the 
study of the Schurness of the matrix Je AT , which turns out to be a very simple 
problem. However, this formulation suffers from several critical drawbacks: 

1. The eigenvalue analysis is not extendable to the case of aperiodic impulses 
since the spectral radius is, in general, not submultiplicative; 

2. To overcome this, Lyapunov approaches can be applied but lead to ro- 
bust Linear Matrix Inequalities (LMIs) with scalar uncertainties at the 
exponential, known to be complex numerically, yet solvable (18|; 

3. The extension to robust stability analysis is also difficult, again due to the 
exponential terms. There is no efficient solution, at this time, to handle 
block matrix uncertainties at the exponential. 

The approach discussed in this paper aims at overcoming the above important 
drawbacks and brings an efficient solution for the characterization of robust- 
dwell times. The main tool used in this paper is an extension of the one de- 
veloped in 20] in the context of sampled-data systems, itself triggered by a 
somewhat different approach discussed in the anterior work [l7| where impul- 
sive systems are considered for the representation of aperiodic sampled-data 
systems 24, 21 , j| . The core of the approach relies on the implicit but equiva- 



lent correspondence between discrete- and continuous-time domains obtained in 
[20| ] . It is shown there that discrete-time stability is equivalent to a certain kind 
of continuous-time stability, provided that the latter is proved using particular 
functionals referred to as looped- Junctionals. One important feature of point- 
wise stability criteria is to allow for non-monotonic continuous-time Lyapunov 
functions (evaluated along the flow a system) since only a pointwise decrease is 
imposed instead of a continuous one. A discrete-time criterion provides then a 
weaker condition for stability than classical continuous-time ones. In the case 
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of impulsive systems, by studying the decrease of the Lyapunov function evalu- 
ated at impulse-times only, non-monotonicity of the continuous-time Lyapunov 
function between impulse-times may be tolerated. 

The interesting point of the method is its wide adaptability to any type of 
systems having discrete-events, or more generally time-marks, which are time- 
instants for which the system admits certain regularity properties. For sampled- 
data and impulsive systems, these time-marks coincide with the control law 
update and impulse-times, respectively. Hence, by looking at a discrete-time 
Lyapunov function whose values at the marker-times are decreasing allows one 
to prove stability of the overall hybrid system. This concept is immediately 
extendable to other type of systems for which such marks may be defined. 

The looped-functional-based approach introduced above leads to LMI-conditions 
that are affine in the dwell-time, devoid of exponential terms and able to con- 
sider the jumps precisely through the non-monotonicity of the Lyapunov func- 
tion. Indeed, expansive jumps on a certain state can be tolerated as long as it 
decreases sufficiently between jumps. Conversely, increasing between jumps is 
also allowed as long as the jumps are contracting. The obtained LMI conditions 
easily allow one to both characterize stability under periodic and aperiodic im- 
pulses. In both cases, the dwell-time Tk = ifc+i — tk plays an important role 
in the stability of the impulsive system. Several stability concepts are consid- 
ered in this paper: stability with minimal dwell-time Tk £ [T m i n , +00) where 
Train is the so-called minimal dwell-time, stability with maximal dwell-time 
Tk £ (0,T max ] where T max is the maximal- dwell time, stability with ranged 
dwell-time Tk £ [T m in,T ma x} and stability with arbitrary pulsing Tk > 0. The 
domains of application of these different dwell-time concepts are summarized in 
Table [TJ It is important to stress that the case (3,3) cannot be handled using 
existing approaches, e.g. d, [l2j], since none of the system matrices are stable, 
while the proposed approach does not require such stability conditions on A 
and J. This emphasizes that the proposed approach better captures the inter- 
nal structure of the system. Thanks to the affine dependence of the obtained 
conditions on the system matrices, dwell-time notions are finally extended to 
the case of uncertain systems, a problem for which there is currently a lack of 
efficient solution techniques. 

It seems important to stress that since the paper has been submitted, several 
improvements have been made. A more advanced explicit looped-functional has 
been proposed in An approach based on an implicit looped-functional has 
also been discussed in Nl. 



Outline: Section [2] introduces a generalization of the result of [20(, at the core 

of the approach. Sections [3] and 2] are devoted to nominal and robust stability 
analysis using a class of looped-functionals fulfilling the conditions stated by 
the main theorem of Section [2] Illustrative examples are included in the related 
sections. 

Notations: For symmetric matrices A,B, A — B -< (^)O means that A — B is 
negative (semi)dcfinitc. The sets of symmetric and positive definite matrices of 
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A 

J ^\ 


U[X(A)} < 


U[X(A)] > 


otherwise 


p(J) < 1 


arbitrary 
minimal DT 


maximal DT 


ranged DT 
maximal DT 


P (J) > 1 


minimal DT 






otherwise 


minimal DT 




ranged DT 



Tabic 1: Application table for the developed results (DT stands for 'Dwell-Time') 



dimension n are denoted by S n and §" + respectively. The spectral radius and 
the spectrum are denoted by p(-) and A(-) respectively. Given a square matrix 
A, we define Sym[A] = A + A T . 

2. Preliminary definitions and results 

2.1. Lifting 

The functional based approach relies on the characterization of system ([1]) 
trajectories in a lifted domain, similar to the one used in sampled-data systems 
theory [26| , with the difference that the involved functions do not have identi- 
cal support. Indeed, we view here the entire state-trajectory as a sequence of 
functions 

{x{t k +t), r e (0,T fc ]} fceN . 

The elements of the sequence have unique continuous extension to [0, Xy defined 

as 

Xfc(r) := x(t k +T), 

Xfc(0) = liWs). ( 3 ) 

We also have the following identities 

Xfe(r) = e AT Xk (0), U] 
Xfc+i(0) = J Xk (T k ). 1 ' 

Hence, in the following, the state-space of the impulsive system will be defined 
as the union set of continuous-functions 



kT mintTmax ] :=U{C([0,T],K«)} 



re [T„ 



with support in a certain range. This varying support is necessary to consider 
the aperiodicity of the system. Note that when T m i n = Tj- aax = T, the usual 
lifting space of periodic sampled-data systems is recovered [26| and degenerates 
to K T := Kpi ,?] = C([0, T], R"). 
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2.2. Stability analysis - Looped- functional-based results 

The definition below introduces the asymptotic stability of impulsive sys- 
tems: 

Definition 2.1. Given an increasing sequence of impulse instants {tk}k£N hav- 
ing no accumulation point, the system (QJ) is globally asymptotically stable if for 
all xq G 1" and all 5 > 0, there exists So > such that 

• ll^oll < So =^ 1 1^(0 1 1 < <5 /or all t > to (stability), 

• ||x(i)|| — > as t — ^ oo (attr activity). 

Alternatively, the asymptotic stability can also be defined in the state-space 
^-[T min ,T max ] where the stability and attractivity definitions become 

• \\ x o\\ < So => sup ||xfe( s )ll < $ f or a M k G 

se[o,T k ] 

• sup ||Xfe( s )ll ^0 as k — > oo. 
s e[o,T fc ] 

It is important to stress that the above stability definition strongly depends 
on the considered impulse sequence and may be very difficult to apply. The 
following weaker definition addresses the case where an entire family of impulse 
instants is considered: 

Definition 2.2. The system is said to be globally asymptotically stable under 
ranged dwell-time if for any increasing sequence of impulse instants {tk}k&i 
having no accumulation point and satisfying tfc+i — tfc G [T m in,T max \, T m i n > 
£ > 0, k G N, the system (QP is globally asymptotically stable. 

The above definitions are clearly not the most general stability notions for 
hybrid systems, but are sufficient for the considered problem. The system being 
linear, properties such as stability are automatically global. Secondly, since 
accumulation points in the sequence of impulse times are excluded, the standard 
asymptotic stability notion is meaningful. It is hence not necessary to define 
hybrid time domains and the positive real axis (or the set of natural numbers) 
can be used to denote time. More general stability notions for hybrid dynamical 
systems, such as pre-asymptotic stability, can be found in [111 ]. 

The following technical definition is necessary before stating the next very 
important result. 

Definition 2.3. A functional f : [0,T max ] x M. [TmzntTmax] x [T mm ,T max ] -> E, 
£ < T m i n < T max < +oo, e > 0, is said to be a looped- functional if the 
following conditions are satisfied 

1. the equality 

f(0,z,T) = f(T,z,T) (5) 
holds for all functions z G C([0, T],R n ) C ^[T min ,T ma;c ] an ^ <M T E 

T m in i Tmax] 7 &nd 
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2. it is differentiable with respect to the first variable with the standard defi- 
nition of the derivative. 

The set of all such functionals is denoted by 2$([T m i n ,T max \). 

The idea for proving stability of ((TJ) is to look at the behavior of a candi- 
date discrete-time Lyapunov function evaluated at the impulse instants, that 
is, we look for a positive definite quadratic form V(x) such that the sequence 
{V (xk(Tk))}ken is monotonically decreasing. This is formalized below through 
a functional existence result: 

Theorem 2.4. Let e < T m i n < T max be three finite positive scalars and V : 
M" — > R + be a quadratic form verifying 

Va;eM™, M1IMI2 < V(x) < (6) 

for some scalars < fj,i < fJ-2- Assume that one of the following equivalent 
statements hold: 

(i) The sequence {V(xk(Tk))}keN is decreasing; that is V(x) is a discrete-time 

Lyapunov function for system x(t k+ i) = e ATk Jx(t k ), T k G [T mm ,T max ]. 

(ii) There exists a looped-functional V G £$([T m i n , T max \) such that the func- 

tional Wk defined for any Ck G K as 

W k (T,Xk,Xk~i) := ^-Ak + V(xk(T))+V(T,Xk,T k )+c k , (7) 
where 

A k = V(xk(0))-V(xk-i(T k -i)), 

has a derivative along the trajectories of system Xk( T ) — Axk{i~), t G 
[0,T fe ] 

^W k (T,Xk,Xk-t) ■= y K * + ^ V ^ Xk ^ + ^V(r,Xk,T k ) (8) 

which is negative definite for all r G (0,7/;), Tk G [T m j n , T max ], k G N. 

Then, the solutions of system (QP with known matrices A and J are asymp- 
totically stable for any sequence of impulse time instants {t k } ke ^ satisfying 
^fc+i tk G [T m j n , T max ] , k G N. 

Proof : Proof of (ii)=>(i): Let k G N, t G (0,T k ] and T k G [T mmi T max ]. 
Assume that (ii) is satisfied. Integrating ((8]) over [0, Tk] yields 

f Tk d 

/ -rWk(T,Xk,Xk-i)dT = W k (T k ,Xk,Xk-i) - Wk{0,Xk,Xk-i) 

= V( Xk (0))-V(xk-i(Tk-i)) 
+V(xk(T k ))-V(x k (0)) 
+V(T k , X k,T k )-V(0,Xk,T k ). 



1 We may also look at the sequence {V(\k (0))}feeN instead. The choice is purely arbitrary. 
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The terms on the last row vanish according to (|5|) and we get 



W k (T k ,Xk,Xk-i)-W k (0,Xk,Xk-i) = V{ Xk {T k ))-V{xk-i{T k -x)). 

Then, the sequence {V {xk{Tk))} km is decreasing over k since ^ is negative 
over[0,T fe ]. 

Proof of (i)=>(ii): Assume that (i) is satisfied. Similarly as in [2(j, intro- 

duce the functional V( T , X k,T k ) = -V( Xk (r)) + [V( Xk (T k )) - V( Xk (0))}. It 

. k 

is immediate to see that V € £$([T m i n , T max \) since we have 

V(T k ,Xk,T k ) = 'V( Xk (Tk)) + (V(xk(Tk)))-V(xk(0)) 

= -V( X k(0)) (9) 
= V(0, Xfc ,T fe ) 

for all Tk G [Tmi n , T max \. Substitution of the proposed functional V into ([8]) 
yields 

^m{r, Xk ,Xk-i) = y [V( X k(T k )) - 7(xfc-i(T fc _i))] 
and is negative by assumption. Equivalence is hence proved. 

Proof of asymptotic stability: It remains to prove that convergence 
of the discrete-time Lyapunov function implies boundedness and asymptotic 
convergence of the continuous-time Lyapunov function, or equivalently of the 
convergence of ||a;(t)||2 to 0. From the discrete-time Lyapunov condition we 
have that ||xfc(^fe)||2 — ^ as k — > oo. Note also that 

V( X k+i(r)) = 

< 
< 

where 77:= sup \(j T e A s Pe As J ] and A is the largest eigenvalue. Note 

that Tmax may be allowed to be unbounded provided that A is Hurwitz since, in 
that case, rj would be bounded. Finally, from the boundedness and convergence 
of V(xk(T k )), we have the boundedness and convergence of V(x(t)) and ||a;(t)||2. 
This completes the proof. 



The interest for considering discrete-time Lyapunov functions lies in the po- 
tential use of non-monotonic continuous-time Lyapunov functions (along the 
trajectories of the system) that are impossible to consider via the usual Lya- 
punov Theorem. Indeed, despite being non-monotonic, Lyapunov functions 
may asymptotically tend to 0, and thus contain information on the asymptotic 
stability of the system. This feature is extremely important in the current 



Xfc+iO) Pxk+i{r) 

Xk{T k ) T J T e ATT Pe AT J Xk (T k ) , . 

ri\\Xk{T k )\\l K ' 

W^VixkiTk)) 
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framework in order to cope with expansive jumps and unstable continuous-time 
dynamics. Using such a discrete-time approach, only the decreasing of the func- 
tion evaluated at instants {tk}keN is important. In Fig. [TJ we may see that 
the two envelopes, generated by the pre-impulses and post-impulses values of 
the continuous-time Lyapunov function, characterize asymptotic stability. The 
functional W of Theorem 12.41 coincides with the lower envelope and a specific 
sequence {ctjfcgs, making the envelope absolutely continuous. Note that the 
lower and upper envelope are equivalent in terms of stability measure since they 
are related through the equality V(xk(fy) — V(Jxk-i(Tk-i))- 

25 
20 
15 

> 

10 
5 

°0 0.5 1 1.5 2 2.5 3 3.5 4 

Time t [sec] 

Figure 1: Continuous-time Lyapunov function V(x) = x T Px (plain) for system l|28| l and 
the discrete-time envelopes (dashed); W coincides with the monotonically decreasing lower 
envelope. 

According to the choice of the sequence {cfej^gNj the functional W can be 
discontinuous. This, however, is not a problem since the functional decrease 
must hold over the intervals (tfc,tfc+i). It is important to keep in mind that 
only the integral of the derivative over (tfc,tfc+i) is important since it coincides 
with a pointwise decrease of the Lyapunov function at impulse instants. 

It will be shown that such an approach is able to characterize stability for 
a larger class of systems than existing methods [H, [l2| , more precisely, systems 
for which neither A nor J are stable. The approach also readily extends to the 
case of uncertain systems for which efficient results are needed. 

3. Nominal stability analysis of linear impulsive systems 

This section provides some nominal stability results for both periodic and 
aperiodic impulses. Therefore, we tacitly assume in this section that the sets A 
and J are reduced to the singletons {A} and {J}, respectively. In the following, 
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we will make extensive uses of the matrix expressions: 



®{P, J) 
S(P,A,J,T) 



= A P + PA, 

= J T PJ - P, (11) 
= J T e ATT Pe AT J - P, 



where ff, 2# and J? stand for 'continuous', 'discrete' and 'impulsive', respec- 
tively. 

3.1. Stability analysis of linear impulsive systems - Periodic impulses case 

The simple case of periodic linear impulsive systems is considered first to 
set up ideas and key results on minimal dwell-time, maximal-dwell-time, arbi- 
trary impulses and ranged dwell-time. In the periodic impulse case, the LTI 
continuous-time impulsive system (fTJ) can be converted into the LTI discrete- 
time system 

x(t k +i) = e AT Jx(t k ), (12) 

the asymptotic stability of which is equivalent to the asymptotic stability of (fTJ). 
This is formalized in the following lemma: 

Lemma 3.1 (Periodic impulses). Assume T > 0, then the following state- 
ments are equivalent: 

a) The LTI continuous-time impulsive system (QP with T -periodic impulses 
is asymptotically stable. 

b) The LTI discrete-time system A12\) is asymptotically stable. 

c) There exists a matrix P = P T >- such that 

S(P,A,J,T) -< 0. (13) 

The main drawback of the above results lies in the presence of exponential 
terms in LMI (fT"3"|) preventing any extension to robust dwell-time characteriza- 
tion, essentially due to the difficulty of considering uncertainties at the expo- 
nential. Note also that the LMI (TTBl is difficult to check when T belongs to a 
(possibly infinite) interval, as this is the case for certain dwell-time notions. 

Below, we show that Theorem l2.4l mav be used to obtain sufficient conditions 
for the asymptotic stability of system ([T]) in the case of periodic impulses. The 
same result turns out to be useful for deriving alternative dwell-time results 
and generalizing them to uncertain systems by resolving the exponential-related 
tractability problem mentioned above. 

Theorem 3.2. The impulsive system (QJ) with T-periodic impulses is asymp- 
totically stable if there exist matrices P, Z E S" + , Q,U G R £ M" xn and 
N € R nx2n such that the LMIs 



*(T) := F +T(F 2 + F 3 ) ~< 

~F -TF 3 N T ~ 
$(T) := Z -< 

* ~T- 



(14) 
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Xk{s) T Zxk{s)ds (16) 
Xk(0) T U X k(0) 



hold with M x = [l 0] , M c = [l -J] , M_ = [0 I], F 3 = M T J T UJM- and 

F = TM x r (A T P + PA)M X - M^QMq 

+M^(J T PJ - P)M_ + Sym[N T M c - M^RM X ] , 
F 2 = Sym[MjA T QM <: + M x A T RM X + M£RAM X ] [ ° } 

+M x r A T ZAM x . 

Moreover, the quadratic form V(x) = x T Px is a discrete-time Lyapunov func- 
tion for system (QJ), that is the LMI ^{P, A, J,T) -< holds, implying then the 
satisfaction of the conditions of Lemma \3.1[ 

Proof : The proof is inspired from [20(. Choosing V(x) — x T Px and 

V(r,x*,T) = ^-fick(r) T [QCk(r)+2R Xk (r)} 
, (T-t) 
T 

where Cfc(r) = Xfc (r) - X fc(0), * fc (r) = A Xfc (r), P g S%,Z g §" + ,Q g 8", C/ e 
S n and i? g R nx ™. The condition © is verified since 

V(0,z,T) = V(T,z,T) = (17) 

for all z g C([0,T],R n ). Thus, according to Theorem COO the functional W k 
defined in (J7J must be considered. Its derivative (|8]) taken along the trajectories 
of the system is bounded from above by the quadratic form 

VV fe < ^Ur) T [Fo + rF 1 + (T-T)F 2 + (T-2T)F 3 }^ k (T) (18) 

where F x = N T Z~ X N and £ fc (i) = col(xfc(r), Xfc-i (7 1 )). The above bound on 
Wk has been obtained using the affine Jensen's bound [3, 0] on the integral 
term 

" Xk (s) T Z Xk (s)(s)ds < ^{t) t {2N t M c +tN t Z- 1 N)^{t) (19) 

defined for some N £ R nx2n . This bound is known to be more precise and 
relevant than the rational Jensen's inequality when applied on integrals with 
uncertain/ varying integration bounds [2|. Hence, the system ([TJ with periodic 
impulses is asymptotically stable if W fe is negative definite over r g (0, T), or 
equivalently, if the LMI 

Fa + tFi + (T - t)F 2 + (T - 2r)F 3 -< 

holds for all r £ [0, T]. Since this LMI is affine in r (hence convex), to check its 
negative definiteness over the entire interval (0, T], it is necessary and sufficient 
to check it at the vertices of the set, that is only over the finite set r g {0,T}. 
A Schur complement on the quadratic term TN J Z _1 N then yields the result. <J) 
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Remark 1 (Note on the affine Jensen's bound). By eliminating the slack- 
matrix N on the RHS of f!9\) using the Finsler's lemma \2ty the usual Jensen's 
bound is recovered, showing then their equivalence. However, as pointed out in 
fl/, they are not equivalent from a computational point of view. Indeed, the 
rational Jensen's bound is nonconvex in the measure of the integration support 
(here t) and must be bounded by its maximal value, i.e. T, resulting then in a 
conservatism increase and a loss of equivalence. On the other hand, the affine 
Jensen's bound is affine in the measure of the integration support (hence convex) 
and does not need to be overbounded in order to make the conditions tractable. 

It is important to stress that while Theorem 12.41 provides a necessary and 
sufficient condition, the above result provides a sufficient one only. By indeed 
choosing the specific functional (fl"6j) . necessity is destroyed. Despite of that, 
it will be shown in the examples that such a functional may yield interesting 
results. To anticipate a little bit, the functional-based formulation will also 
be shown in Section [5] to be suitable for robust stability analysis. All these 
characteristics emphasize the relevance of the proposed approach inheriting part 
of the efficiency of the usual discrete-time Lyapunov approach. 

We digress here a little bit to address the case of stability under small periods, 
i.e. the case T — > 0, but excluding any Zeno behavior. 

Lemma 3.3. When T — >• 0, then the LMI conditions \1J$ tend to the condition 
@(P, J) -< 0, equivalent to the Schurness of J. This is consistent with the fact 
that p{Je AT ) -> p(J) asT -> 0. 

Proof : When T -» 0, then we have *(T) -)■ F | T=0 and $(T) -)• F | T=0 where 
Fo\ T=Q is given by 

-M C T QM C - M^RM X - MjR T M c 

+N T M C + M^N + M T (J T PJ - P)M- [ ' 

and should be negative definite. By virtue of the Finsler's Lemma [2^, the 
matrix TV can be eliminated and we obtain the equivalent LMI 

MlM^(J T PJ - P)M-M± -< (21) 
where M± — [j T l\ T . Evaluation of the expression yields the result. () 

It is important to stress that T — > does not necessarily means that an 
accumulation point exists. The limit in the small period allows one to study the 
stability for very small periods that are still different from 0. Let us consider 
for instance the unbounded sequence of impulse times given by t^ = log(fc + 1), 
k 6 N. The dwell-times are then given by Tf. := tk+i — tk — log(l + l/(fc + 1)) 
and tends to as k goes to infinity. In this case, stability of the system must 
be ensured for arbitrarily small positive dwell-times. 

Thus, from Lemma 13.31 and by virtue of the continuity of eigenvalues with 
respect to matrix parameters, the condition p(J) < 1 implies p(e AT J) < 1 in a 
sufficiently small positive neighborhood of T = 0. This result will turn out to 
be very useful in Section 13.41 on the ranged dwell-time. 
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3.2. Minimal dw ell-time characterization 

Several results characterizing stability with minimal dwell-time and arbitrary 
impulse sequences are provided below. The two results considering the minimal- 
dwell-time are obtained using similar arguments as in [Toj . and using Theorem 
12.41 respectively. The case of arbitrary impulse sequences is a particular case of 
the first result. 

Lemma 3.4 (Minimal Dwell-Time). Assume that for some given T > 0, 
there exists a matrix P — P T y such that the LMIs 

tf(P,A)^0 (22) 

and 

J{P,A,J,T) <Q (23) 

hold. Then, for any impulse sequence {tk]k&i satisfying t k+ i — t k > T, the 
system (Qp is asymptotically stable. 

Proof : The goal is to show that the conditions ([22]) and (|23]l implies that 
J^(P, A, J, s) -< for all s >T. To this aim, we use the fact that the eigenvalues 
of J T e A s Pe As J are nonincreasing as s > increases. To show this, let us 
consider the Lyapunov function evaluated at t k + s, s > given by 

V(x(t k + s)) = x(t k ) T J T e ATs Pe As Jx{t k ) (24) 

whose derivative with respect to s is given by 

^-V(x(t k + a)) = x(t k ) T J T e ATs tf(P, A)e As Jx(t k ). (25) 
as 

By virtue of the condition (f2"2"j) . the derivative is negative scmidcfinitc for all s > 
and hence the eigenvalues of J T e A s Pe As J are nonincreasing as s increases. 
We then have the inequality 

J T e ATs Pe As J d J T e ATT Pe AT J 
for all s >T. Using now the condition (|2"3")l this implies that we have 

■y(P,A, J,s) ^ J?{P, A, J,T) -< (26) 
for all s >T. The proof is complete. ()> 

Similarly as in the periodic case, we derive below an affine counterpart of 
the above result, which will be useful for the minimal dwell-time analysis of 
uncertain systems. 

Theorem 3.5 (Minimal Dwell-Time). Assume that for some given T > 0. 
there exist matrices P, Z e §I^ + , U,Q € S n , N G R nx2n and R G M" xn such 
that the LMIs *(T) -< 0, $(T) -< and tf(P,A) -< hold. 

Then, for any impulse sequence {t k } ke jq satisfying t k +i —t k >T, the system 
(QP is asymptotically stable and Lemma \3.J\ holds with the same matrix P. 
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Proof : The proof follows from Lemma 13.41 and Theorem 13.21 (} 

Note that when deriving Lemma l3.4[ there is no restriction in letting T Q, 
leading then to a possible characterization of the stability for arbitrarily small 
dwell-times. Lemma 15731 savs that if there exists P G §" + such that J, P) -< 
holds, then asymptotic stability is guaranteed for sufficiently small dwell-times. 
The following corollary is a consequence of this fact and the maximal dwell-time 
result of Lemma 13.41 

Corollary 3.6 (Arbitrary impulses sequence). The system {J) is asymp- 
totically stable for arbitrary impulse sequence {tfejfeeN verifying tk+\ — tk > e, 
for any < e < T and for all k G N, if there exists a matrix P G S" + such that 
C(P, A) -< and 3>{P, J) -< hold. 

Proof : It is sufficient to let T — > in the condition (|2"3")) of Lemma 13.41 Q 

Note that the same condition can be found in [27], Theorem 4.1.1] or in [l2l 
Theorem 2] albeit expressed in different ways. 

3.3. Maximal dwell-time characterization and arbitrary impulses 

Results on stability with maximal dwell-time [12( are discussed in this sec- 
tion. 

Lemma 3.7 (Maximal Dwell-Time). Assume that for some given T > 0, 
there exists a matrix P G §" + such that the LMIsC(P, A) >~ and J{P, A, J, T) -< 
hold. 

Then, for any impulse sequence {£fc}fc<=N satisfying e < t^+i — tk < T, for 
any < e < T , the system (Qp is asymptotically stable. 

Proof : The proof is similar to the one of Lemma 13.41 {> 

Theorem 3.8 (Maximal Dwell-Time). Assume that for some given T > 0, 
there exist matrices P, Z G y 0, U, Q G S", N G R nx2n and R G W ixn 
such that the LMIs *(T) -< 0, $(T) -< and tf(P, A) y hold. 

Then, for any impulse sequence {tk}kefi satisfying e < tk+i — tk < T, for 
any < e < T , the system f7|) is asymptotically stable and Lemma \3. 7| holds 
with the same matrix P. 

Proof : The proof follows from Lemma 13.71 and Theorem 13.21 ^ 

3.4- Stability over arbitrary intervals - Ranged dwell-time 

The above results can only be applied when A is Hurwitz or anti-Hurwitz. To 
overcome this limitation, the following results considering dwell-times belonging 
to general bounded intervals of the form [T m i n ,T max } are provided. This leads 
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to the notion of ranged dwell-time which combines the concepts of minimal and 
maximal dwell-time together. 

It is immediate to derive the following preliminary result on ranged dwell- 
time: 

Lemma 3.9 (Ranged dwell-time). Assume there exists a matrix P G §™ + 
such that 

S(P,A, J, 9) -< (27) 

for all 9 G [T min ], e < T mm < T max < +oo, for some e > 0. 

Then, for any impulse sequence {t k } keti satisfying tk+i - t k G [T min ,T max ], 
the system (QP is asymptotically stable. 

It is clear that the semi-infinite dimensional LMI (|27p is not easy to check due to 
the presence of the uncertain parameter 9 G [T m i n ,T max \. The use of Theorem 
12.41 allows one to overcome this difficulty 

Theorem 3.10. The impulsive system f7]] with T k G [T m i ni T max ], e < T m i n < 
T m ax < oo, e > 0, is asymptotically stable if there exist matrices P, Z G 
Q, U G §™, R G R nxn and N G R nx2n such that (T) -< and $(T) -<! /iota 
/or cz/£ T G {T m 2n , T 7nax } . 

Moreover, in such a case, the inequality J'iP, A, J, 9) -< ZioZrfs /or 
G [T m i n ,T max ] and Lemma \3.9\ is verified with the same matrix P. 

Proof : Let us consider the LMIs of Theorem 13.21 allowing the check the sta- 
bility of the system for a constant dwell-time T > 0. In order to prove that the 
impulsive system is asymptotically stable for a varying dwell-time belonging to 
[T m in,T max ], the LMI must be simply checked over the interval 9 G [T m i n , T max \. 
Since the LMIs \&(T) and $(T) are convex in T, it is then enough to check it at 
the vertices of the interval, and hence for all values in the finite set {T m i ni T max }. 
This concludes the proof. (} 

It seems important to point out that, unlike the results on minimal and 
maximal dwell-time requiring stability or anti-stability of the matrix A, the 
above result is applicable to any linear impulsive system. It can notably be 
applied to systems for which neither A nor J is stable, which is very interesting 
since existing methods cannot handle such a situation. This is to put in contrast, 
for instance, with |12l Theorem 1], transferred to a linear setting, for which it 
is necessary that at least one of the matrices be stable. Another feature is the 
reduction to a finite-dimensional LMI problem while the condition of Lemma 
13.91 is semi-infinite dimensional. 

Notably, the above result may be used to provide an alternative maximal 
dwell-time result in which the anti-stability constraint on the matrix A is re- 
laxed: 

Corollary 3.11 (Alternative maximal dwell-time result). Assume that for 
some given T > 0, there exist matrices P, Z G S" + y 0, U,Q £ S n , N G R nx2n 
and R G K" xn such that the LMIs *(T) -< 0, $(T) -< and 9{P, J) -< hold. 
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Then, for any impulse sequence {tk}ken satisfying e < tk+i — tk<T,e> 0, 
the system (QJ) is asymptotically stable and Lemma \3.7\ holds with the same 
matrix P. 

Proof : Using the ranged dwell-time result (Theorem 13. 10[) and letting T — > 0, 
we obtain the limit LMI J T PJ — P ~< 0, as shown in Lemma I3~3"l This concludes 
the proof. (} 



3. 5. Examples 

It is illustrated below, through academic examples, that the proposed alter- 
native dwell-time formulations may lead to interesting results in terms of accu- 
racy. It is important to stress that the conservatism is only due to the choice for 
the functional since Theorem 12.41 provides equivalent statements. This increase 
of conservatism is the price to pay to obtain tractable and efficient tools for the 
analysis of uncertain linear impulsive systems. 

Example 1 (Maximal dwell-time). Let us consider system $7$) with matri- 
ces 

' 1 3 



.4 



-1 2 



J 



0.5/ 2 



(28) 



Since A has 2 unstable eigenvalues at 1.5 ± 1.6583j, then overall system stability 
requires an asymptotically stable J , which is the case here. Hence if the pulses 
do not occur frequently enough, stability is not achievable. 

Periodic impulses case: In the periodic case, the spectral radius condition 
can be used and is equivalent to the feasibility of J? (P, A, J, T) -< for some P G 
S" + . The spectral radius condition gives the maximal impulse period T max = 
21og(2)/3 G (0.4620,0.4621). Theorem HOI is applied together with a bisection 
approach to find the maximal constant T > that preserves stability and yields 
the lower bound T max = 0.4471. This shows that the provided approach is able 
to estimate quite well the maximal pulse period for this example. 

Choosing for instance T = 0.3, we obtain the states and continuous-time 
Lyapunov function trajectories for system i28\) depicted in Fig. [H and [7J 

Aperiodic impulses case: Since the matrix A is anti-Hurwitz, then this 
system may fulfill conditions of the maximal dwell-times results. When Lemma 
\3. 7[ based on the discrete-time Lyapunov condition, is considered, the lower 
bound on the maximal dwell-time T^ aa , = 0.4620 is found. This shows that 
Lemma \3.7\ is able to exactly characterize the maximal dwell-time since the com- 
puted value is identical to the maximal admissible impulse period (periodic case). 
The system 1128]) is hence asymptotically stable for all Tk G (0, T max ] . Theorem 
\3.8\ is now considered and yields the lower bound on the maximal dwell-time 
Tmax = 0.4471. Hence, Theorem \3.8\ vrovides. for this example, a good approx- 
imation of Lemma \3. 



Example 2 (Minimal dwell-time). Now consider system {IJ) with matrices 



A 



-1 
1 



J 



2 1 
1 3 



(29) 
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Time t [seel 



Figure 2: States trajectory of impulsive system \28) with impulse period T = 0.3. 



Since the matrix A is Hurwitz and J is anti-Schur, then, if the impulses are 
too frequent, stability is lost. It is hence expected to find a sufficiently small 
dwell-time for which stability is lost. 

Periodic impulses case: The spectral radius condition yields T TO , n £ 
[1.1405,1.1406] while Theorem \3.S\ gives the upper-bound T^ in = 1.2323 on the 
minimal impulse period. 

Aperiodic impulses case: Since A is Hurwitz, then the system may ful- 
fill conditions of the minimal dwell-time results. Lemma \3.4\ yields the upper 
bound on the minimal dwell-time 1.1405 which is almost identical to the min- 
imal period of the periodic case. This emphasizes the exactness of Lemma \3.4\ 
in the estimation of the minimal dwell-time. System i29\) is hence guaranteed 
to be asymptotically stable for all Tk G [T^ in , +oo). In contrast, Theorem \3.5\ 
predicts the upper bound T^ in = 1.2323, showing the quite good accuracy of the 
proposed approach. 

Example 3. Consider now system (QP with matrices 



A = 



-1 




0.1 
1.2 



J = 



1.2 






0.5 



(30) 



Note that, the continuous-time dynamics of the first state is stable while the 
second is unstable. Conversely, the matrix J has a stable eigenvalue for the 
second state and an unstable one for the first state. It is hence expected that 
the range of admissible dwell-times is a connected interval excluding and +oo. 
It is important to stress that such a system cannot be analyzed using existing 
methods, such as the ones in 0, [71/ . since neither A nor J is stable. 

Periodic impulses case: An eigenvalue analysis gives the admissible range 
[0.1824,0.5776] of dwell-times. When Theorem \ 3. 6 A is used with two bisection 



routines, the following bounds are obtained T T ^ i 



0.1824, Ti 



0.5760. 
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This illustrates that the approach is able to characterize quite well the stability 
of this system. 

To emphasize that this case cannot be handled with existing methods, let 
T = 0.3 and we find P = diag(2.3622, 1.4752) along with fH]: 



,A) d 
@{P, J) < 



-cP 



DP 



(31) 



where c = —2.4036 and d = —0.3646. Note also that no P y satisfies one of 
the conditions &31\) with c > or d > 0. Since both c and d are negative for all 
P y 0, the method of fl$q is inconclusive. This demonstrates the potential of 
the proposed functional-based approach. 

Aperiodic impulses case: Theorem VJ. 1 0\ yields the interval [0.1907, 0.5063], 
which is included in the interval obtained for the case of periodic impulses. 



This example demonstrates that the proposed approach is able to capture 
the internal structure of this system better than other approaches, and take the 
advantage of this feature to improve accuracy. The main drawback of the existing 
methods based on a Lyapunov conditions and a-stability of the form i31\) is that 
they consider the worst case eigenvalue (covering of the system) and lose/ignore 
the internal structure of the system, notably eigenvectors. 



Example 4. Let us consider now the sampled-data control system 

x(t) = Ax(t) + Bu(t), 

u(t) = Kx(t k ), t e [t k ,t k+ i) 

Reformulating this system as an impulsive system, we obtain 



(32) 



m 

z + {t k ) 



A B 



In 0' 

K 



z(t) 
z(t k ) 



(33) 



where z(t) — col(x(t), u(t)). Note that neither A nor J is a stable matrix, hence 
the developed maximal and minimal dwell-time results cannot be applied. The 
ranged dwell-time result however applies here. It is well-known that when the 
corresponding continuous-time system is asymptotically stable, i.e. A + BK 
Hurwitz, then the sampled-data system is stable in a sufficiently small positive 
neighborhood of the 'zero sampling-period' fl^, [jfj /. Unfortunately, this re- 
sult cannot be used in the ranged dwell-time result (Theorem \ 3.10\) since this 
condition cannot be obtained using the current approach that requires J to be 
Schur. 

Let us consider now a sampled-data system with matrices \v\ , Q/ 



A = 



1 

-0.1 



B = 





0.1 



K = 



[3.75 11.5] . (34) 
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Applying Theorem \3.2\ we find the lower bound on the maximal constant 
sampling period T^ax = 1-7239 similarly as in Q/. Comparatively, the 
spectral radius condition yields a maximal value of 1.7294. Focusing now on the 
aperiodic case, we apply Theorem \3.10\ where we set T m i n = 10 and we find 
the lower hound Tz. n „ = 1.7239. 



4. Quadratic stability analysis of aperiodic uncertain linear impulsive 

systems 

Extensions of the previous results to the case of uncertain systems are dis- 
cussed in this section. We hence now consider the case where the matrices A 
and J belong to some distinct polytopes A and J , as defined in ([2]). Due to the 
affine structure of the conditions, the extension to uncertain systems is imme- 
diate. Despite of this ease, most of the efficiency of the approach is preserved 
as demonstrated in what follows. 



4--1- Main results 

The cases of periodic impulses and minimal, maximal and ranged dwcll-times 
for uncertain systems are discussed. 

Theorem 4.1 (Periodic impulses). The uncertain linear impulsive system 
{IP-{IP with periodic impulses is asymptotically stable if there exist matrices 



P,Zj e Qj,Uj e S n , Rj e 

that the LMIs 



and Nj e 



for j = 1, 



G%+T{G 



Gi 



MlAfZ, 
-Zi/T 



- TGi 



Ni 



hold for all i = 1, . . . , Na and j = 1, 



-Z^T 
. , Nj where 



-< 



-< 



0. 





, Nj such 



(35) 



/~i0 



Gh 



T ■ Sym[MjAfPM x ] - M ' Q,.\l, 
+MT(jTpjj - p)M_ + Sym^/M,- - MjRjM x ], 
Sym [Mj Aj Qj Mj + MjAfRjM x + Mj RjAiM x \, 



G] = M^jfUjJjM- 

Moreover, the quadratic form V(x) = x T Px is a quadratic Lyapunov function 
for the uncertain system {JJ), that is, we have ^{P, A, J, T) -< for all (A, J) E 
Ax J. 

Proof : The goal is to show that the feasibility of the LMIs ^y(T) -< and 
$ij(T) -< 0, for alH = 1, . . . , Na and all j = 1, . . . , Nj implies the feasibility of 
J(P, A, J, T) -< for all (A, J) eAxJ. 

N A 

To show this, multiply first and by K i where k = col[/c,] belongs to 

i—l 

N A N a 

the unit iV^-simplex. The sums ^l^ij and K l^ij are then both negative 

i=i i=i 
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definite since the Vf^-'s and 3>ij's are individually negative definite. After a Schur 
complement, we obtain an LMI of the same form as in Theorem 13.21 implying 
then that for each j = 1, . . . , Nj we have 

J{P,A,Jj,T) -< 

for all A £ A. Note that, by virtue of the Schur complement formula, the above 
LMI is equivalent to 



-P jje AT P 
* -P 



< (36) 



Nj 

which is linear in Jj. Multiplying the above LMI by n 2 where k 2 — co1[k 2 ] 

belongs to the unit TVj-simplex and summing over j = 1 , . . . , Nj yields the 
condition 

J?{P, A, J, T) < 0, for all (A, J) g A x J (37) 
by using again the Schur complement formula. This concludes the proof. 

Following the same reasoning as in the previous section, the above theorem 
can be used to easily derive dwell-time results for uncertain systems. 

Theorem 4.2 (Robust minimal dwell-time). Assume that for some given 
T > 0, there exist matrices P,Zj g §™ + y 0, Uj,Qj g Nj g R nx2n and 
Rj g W ixn such that *ij(T) -< 0,. $y(T) -<! and tf{P,Ai) -< hold for all 
i = 1 , . . . , Na and all j = 1 , . . . , Nj . 

Then, for any impulse sequence {tk}ken satisfying tk+i — tk>T, the uncer- 
tain linear impulsive system is asymptotically stable. That is we have 
J(P, A, J, 0) -< for all {A, J) g A x J arid aZZ 6» g [T, +oo). 

Theorem 4.3 (Robust maximal dwell-time). Assume that for some given 
T > 0, there exist matrices P, Zj g §™ + >- ; Uj,Qj £ S n , Nj £ R" x2n and 
Rj g R" x ™ suc/i i/iai * y (T) -< 0, $ij(T) -< and ^{P,Ai) >- Aoid /or a// 
i = 1 , . . . , Na and all j = 1, . . . , Nj . 

Then, for any impulse sequence {tk}ken satisfying e < tk + i — tk < T, e > 0, 
the uncertain linear impulsive system is asymptotically stable. That is 

we have J{P, A, J, 9) -< for all (A, J) £ A x J and all 9 g (0, T], 

Theorem 4.4 (Robust ranged dwell-time). Assume that for some e < T m i n < 
T m ax < +oo, £ > 0, there exist matrices P,Zj g Qj,Uj £ §™, Rj £ 

R nxn and Nj £ R" x2n such that ^ij{T) -< and $«(T) -< /io/d /or a// 
7 1 g {T mm , T max } ; all i — 1, . . . , -/Va and all j = 1, . . . , AT/. 

T/ien, /or any impulse sequence {tk}keN satisfying tk+i - tk £ [T mm ,T ma x], 
the uncertain linear impulsive system fH)-@j is asymptotically stable. That is 
we have J(P, A, J, 9) -< for all (A, J) £ A x J and all 9 £ [0, T] . 
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4-2. Examples 

Example 5 (Robust Maximal Dwell-Time). Let us consider an uncertain 
version of the system treated in Example^ Assume that the matrix J is known 
and that the matrix A belongs to 



A := co 



1 3 
-1 2 



2 2 
6 



Periodic impulses case: An eigenvalue analysis performed over the set 
of all possible systems locates the maximal impulse period inside the interval 
[0.1155,0.1156]. By gridding the LMI with a very thin grid S(P, A, J,T) <0, we 
obtain the lower bound 0.1155 on the maximal admissible period, at the expense 
of a very high computational cost, gridding imprecision and large computation 
time. In contrast, Theorem \4-l\ yields the lower bound on the maximal period 
Tmax = 0.1148 at a much lower computational cost and infinite precision (since 
the entire interval of dwell-times is considered rather than a finite subset of it), 
but with a slight conservatism. 

Aperiodic impulses case: By gridding the conditions of Lemma \3.7\ the 
lower bound 0.1155 on the maximal dwell-time is obtained. When Theorem \4-4\ 
is used, we obtain T^ nax = 0.1148. This demonstrates the interest of functionals 
for analyzing the stability of uncertain linear impulsive systems. 

Example 6 (Robust Ranged Dwell-Time). Example\^is revisited here with 
the difference that A is known while J G J where 



J := co 



1.3 
0.25 



1.1 







0.5 



An eigenvalue analysis done over the entire set of possible systems yields the 
following interval of admissible impulse periods [0.2624,0.5776]. In contrast, 
Theorem \4 .l\ yields the upper and lower bounds on the admissible impulse periods 
Tmin — 0.2625 and T,^ ax = 0.5761. These bounds immediately generalize to the 
aperiodic case. 



5. Conclusion 

A loopcd-functional-based approach has been proposed for the analysis of un- 
certain linear impulsive systems. The main advantage of the proposed approach 
lies in the expression of a discrete-time stability criterion using a continuous- 
time approach. Such a framework allows one to easily consider non-monotonic 
continuous-time Lyapunov functions and is able to analyze a wider class of 
systems than existing approaches. As a byproduct, this approach leads to an 
exponential-free formulation of the stability conditions, allowing to obtain ro- 
bustness result very easily. Several examples illustrate the potential of the 
approach. 

The proposed methodology can be applied to many other types of continuous- 
time and hybrid systems for which a recursive discrete-time stability criterion 
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can be applied. Furthermore, since the current formulation involves continuous- 
time data, the framework can be readily extended to nonlinear systems, higher 
order Lyapunov functions, the use of T-dependent Lyapunov functions, etc. 
These problems are way beyond the scope of this paper and are left for future 
works. 
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